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Abstract 

We consider a piezothermoelastic panel occupied by a material of 
hexagonal crystal class. We study the response when the boundary 
conditions vary very slowly with time and one of the bounding faces is 
subject to thermal exposure. We show that in some cases the temper- 
ature on the other bounding face can be controlled by the difference 
of electric potential between the faces. 

Keywords: Piezothermoelasticity, Electromagnetic effects, Thermal ef- 
fects, Plates, Piezothermoelastic plate. 

1 Introduction 

1.1 Premise 

Piezoelectricity is the property of generating an electric field (mechanical 
stress) in response to an applied mechanical stress (electric field); pyroelec- 
tricity is the property of generating an electric field (temperature change) by 
a temperature change (electric field). 

For a piezothermoelastic body a natural problem, useful for practical 
applications, is to study the electromechanical effects due to a prescribed 
temperature on part of the boundary. 

The present paper is a starting study of boundary control problems in 
plate-like bodies exhibiting pyroelectricity. This in order to theoretically 
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establish whether and under which boundary conditions a given boundary 
temperature and/or electric potential can be usefully employed to obtain 
some type of control, e.g. on temperature and/or electric potential, in a 
given material surface. This theoric study may be useful with regard to real 
panels subject to sun exposure or lying in contact with an external heat 
source, in order e.g. to passively exploit this boundary condition. 

Here we define first a general boundary control problem and then study 
the particular problem of a plate occupied by a material exhibiting piezo- 
electric and pyroelectric properties. 

1.2 Static Boundary Value Problem for a Piezother- 
moelastic body 

Here we adopt the summation convention and comma notation for partial 
derivatives, so that a l bi = Y2i a% ^i an d f,% = df/dxi. 

Consider a piezothermoelastic body B that, in the reference configuration, 
occupies a region V with boundary surface S. The deformation of the body 
is described by 

Vi = yi (X A ) = 8fX B + Ui(X A ) , 

where denote the spatial coordinates and X L the reference coordinates of 
material points with respect to the same Cartesian coordinate system. 

The Piola-Kirchhoff stress tensor, electric displacement vector, and heat 
flux vector are respectively given by constitutive functions 

tu = tu(9, Et, Eab,Xa) , (1) 

Dl = D L (9, Ei, Eab,Xa) , Ql = Ql{9, Ei, Eab,X a ) , (2) 

where 9, 0, Ei = —<f> t i, Eab = (yi,AUi,B — Sab)/2 are the absolute temper- 
ature, electric potential, electric field and strain tensor, respectively. 

Balance law of linear momentum, Maxwell's equation, and balance law 
of conservation of energy, respectively lead to the equilibrium relations 

tu,L + Pofl = , D L)L = p e , g i i= p o7 , (3) 

where p D is the mass-density in the reference configuration, j) is the body 
force per unit mass, p e is the body free charge density, and 7 is the body 
heat source per unit mass. 

To describe the corresponding boundary conditions, three partitions (Sn, S i2 ), 
i = 1, 2, 3, of the boundary surface S = dB can be assigned. For mechanical 
boundary conditions, displacement u and traction t per unit undeformed 
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area are prescribed, respectively, on Sn and Si 2 ; for electric boundary con- 
ditions, electric potential and surface-free charge D per unit undeformed 
area are prescribed, respectively, on S21 and S22', while for thermic bound- 
ary conditions, temperature T and normal heat flux q per unit undeformed 
area are prescribed, respectivelyon and S32. Hence, we can write 

Ui = Ui on Su , tuNi = ti on S12 ('mechanical') , (4) 

= on S21 , DlNl = D on S22 , ('electric') (5) 

T = T on 5*31 , q L N L = q on S 32 ('thermic') , (6) 

s ll us i2 = s, S a nS i2 = ® (z = l,2,3), (7) 

where N = (Nl) is the unit exterior normal on S and T = — O is the 
incremental temperature with respect to the temperature 8 in the reference 
state. 

The boundary value problem is then stated as: to find the solution (0, T, u) 
in B to the constitutive relations (TJP, (HP and field equilibrium equations ^ 
which satisfies the boundary conditions for given u^, t{, 0, D, T, q. 

Of course, existence and uniqueness of the solution must be separately 
examined, but here are assumed. 

1.3 Boundary Control Problem 

Let S be a regular surface contained in B := B U dB (possibly S C dB), 
with oriented normal unit vector N, let C : S Q —* M be a scalar smooth 
function, and let x be any one of the ten quantities 

0, 6, D L N L , q L N L , u h t Li N L (1 = 1, 2, 3) . 
The boundary control problem for \ on S with goal C, (abbreviated to 
BCP x\s = C)j is the following: For arbitrary n (< 10) boundary conditions 
of the type choose the remaining 10 — n boundary conditions in Q)- 

(EP such that the solution (0, T, u) to the corresponding boundary value 
problem yields x\s a = C- 
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1.4 The Boundary Control Problems solved here 



Boundary Control Problems for panel V 



BCP 7.1.3 


BCP 77.1.3 


BCP 7.3.3 


BCP 77.3.3 


X\ thickness d. 
x 3 polariz. d. 


Xi thickness d. 
x 3 polariz. d. 


x 3 thickness d. 
x 3 polariz. d. 


X3 thickness d. 
x 3 polariz. d. 


at xi — h 
T = T 

= 

tN = t 


at x\ — h 
T = T 

= 

tN = t 


at x 3 = h 
T = T 

= 

tN = t 


at x 3 = h 
T = T 

= 

tN = t 


at x\ = —h 
u = u 
D ■ N = D 
q • N = q 


at x\ = —h 
u = u 

= 02 

q- N = q 


at £3 = —h 
u = u 
D-N = U 
q ■ N = q 


at x 3 = —h 
u = u 

= 02 

q- N = q 


<f>(xi) controlled 
T(xi) controlled 
—h < x\ < h 


T(xi) controlled 
—h < x\ < h 


0(2:3) controlled 
T(xs) controlled 
—h < x 3 < h 


T(x 3 ) controlled 
—h < x 3 < h 



In the present paper, we solve some boundary control BCP, for certain 
linear piezothermoelastic bodies, that occupy a plate V infinite in extent and 
bounded by two parallel planes. The plate has a natural equilibrium state, 
i.e., with no initial field, and is occupied by a heat-conducting piezoelectric 
material with the symmetry of the hexagonal crystal class C Qu = 6mm, so 
that ferroelectric ceramics are included. We assume that V is subject to 
a constant temperature on the upper face that in effect may vary slowly 
with time. On the lower face, the displacement is prescribed, as when, for 
example, V is welded to a fixed flat body. 

We study processes which are homogeneous on each plane parallel to 
the boundary planes, that is, they depend only on the thickness coordinate, 
and, moreover, vary very slowly with time. The precise equilibrium boundary 
value problems studied are summarized in the table and are completely solved 
once their exact solutions are determined. 

In these problems, we take n = 9 and S a any fixed plane parallel to the 
plane boundaries of the plate. In BCP.s 7.1.3 and 7.3.3 either x = T or 
X = 0, i-e., either temperature or electric potential can be controlled on S Q . 

In BCP.s 77.1.3 and 77.3.3 we have x — T, i.e. whatever temperature 
is prescribed at the upper face, the temperature can be controlled on S a by 
the electric potential difference between the two bounding planes. 
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2 Linear Piezothermoelasticity 



2.1 Linear constitutive equations 

The linear constitutive equations are specified below in terms of the constitu- 
tive coefficients: Ckiij = elastic moduli; e^i = piezoelectric moduli; (3ki = 
thermal stress moduli; K kl = dielectric susceptibility; q k = pyroelectric 
polarizability; e kl = permittivity moduli; km = Fourier coefficients; 7 = 
heat capacity; r\ = entropy at the natural state; p Q = mass-density at 
the natural state. These coefficients, each assumed to be constant, satisfy 
the following symmetry conditions: 

C-klij Cijkl Qfeijr Cklji j &kij &kji j (8) 

Pij — (3ji i K-M — K>lk j K>kl = K lk ■ (9) 

With respect to a natural reference state, i.e., a state free from mechanical 
and/or electric fields, and with constant temperature 9 D , we assume the fol- 
lowing standard constitutive equations [I]- [SI P-122], [0]), respectively for the 
stress tensor, electric displacement vector and heat flux vector: 

tki = CkiijU { j — eikiEi — PkiT , (10) 
D k = ekijUij + e ki Ei + u k T } q k = - K kl T t - fcfj-Ej, (11) 

where Ei = —<ft t j , T = 9 — o is the incremental absolute temperature 
with respect to the absolute temperature 9 in the natural reference state. 

2.2 Field equations of equilibrium 

The linearized field equations of equilibrium, obtained by substituting the 
constitutive equations (jTU]) - (TTTl) in the balance laws taking f k = p e = 
7 = are given by 

CkiijU iJk + e i:i i(f) tij - (3 k iT k = (I = 1, 2, 3) , (12) 



CkjiU j>ik - e kj (f) :jk + ti k T ,k = 0, -K k j T ,jk + K>f k <f> tjk = 0. (13) 
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2.3 Use of compressed notation and matrix arrays 

As is well known, the matrix notation consists of replacing ij or kl by 
p or q, where i, j, k, I take the values 1, 2, 3 and p, q take the values 
1, 2, 3 ,4, 5, 6 according to the following relations: 



ij or kl 


11 


22 


33 


23 or 32 


31 or 13 


12 or 21 


p, q 


1 


2 


3 


4 


5 


6 



By virtue of the above identification, the constitutive equations become 

tp CpqSq €{p E{ ftpT , (14) 

A = e iq S q + e ik Ei + &.T, $ = —k u T x - k% E t , (15) 

where Sy = | (u^j + <% = S p when i — j — p — 1, 2, 3, 2SV,- = 

when i j, p = 4, 5, 6, and = e ip . 

3 Hexagonal materials 

3.1 Constitutive equations for ferroelectric ceramics 

The polarized ferroelectric ceramics have the symmetry of a hexagonal crystal 
in class C§ v = 6mm. Choosing x% in the polarization direction, assuming 
j3ij = for i ^ j and putting 

(k = 1,2,3), [AJ=[/3i, fa, Ps, 0, 0, 0], (16) 

the constitutive equations (JTU1) and (JTT}) i become (cf. e.g. [6], [3 p. 58], and 

H!) 

t\ = CnM^i + Ci 2 U 2 ,2 + Ci3«3,3 + e 3 l 0,3 - AT , 
t 2 = Ci 2 Ml,l + CiiU 2) 2 + Cl 3 M3,3 + e 3 l 0,3 ~ ftT , 

*3 = c 13 (u M +m 2 , 2 ) + c 33 M3 i3 + e 33 i3 , (17) 

t 4 = C44 («3, 2 + U 2 , 3) + ei5 0, 2 , 
t 5 = C 44 (m3,1 +«1, 3 ) +615 0,1 , 

*6 = c 66 (« lj2 + 1*2,1) , with cea = (cu - c X2 )/2 , 
-Di = ei 5 (u 3)1 + M 1( 3) -£n 0,i + a^T, 

-D 2 = ei 5 (« 3 ,2 + «2 )3 ) - £11 0,2 + u 2 T , (18) 
-D 3 = e 3 i (u M + u 2>2 ) + e 33 n 3i 3 - £ 33 0, 3 + & 3 T . 
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3.2 Field equations of equilibrium 

When the constitutive relations ( JT71) . ( fTHl) are substituted in the equilibrium 
field equations (TL2l)-(jT3l) we have 

- PlT t i + C n U 1A1 + (C12 + C 66 )m 212 + (Ci 3 + C44)m 3i i3 + 

+ c m u 122 + c 4 4U lj3 3 + (e 3 i + 615)0,13 = (19) 
- P2T2 + e 66 w 211 + (e 66 + c 12 )u l l2 + c n u 222 + (613 + 644)^3 23 + 

+ 644^,33 + (e 3 i + eis)^, 23 = (20) 

- P3T3 + c 44 %, 11 + (e 44 + ci 3 )m 1)31 + 644%, 22 + (644 + 613)^3,23 

+ C 33%,33 + e 150,ll + ei 5 0, 22 + 633^,33 = (21) 

e i5%,n + (ei 5 + e 3 i)« 1)13 + ei 5 u 3i22 + (e 15 + e 3 i)u 2 , 32 + 

e 33%,33 + ^l T J - £ U<f>,U - £110,22 - ^330,33 = (22) 

-><,J.,j + 4<f> tij = 0. (23) 

Particular forms of the solution to these equations are discussed in the next 
section, while general expressions are derived in the Appendix. 

4 Quasi-statics 

A principal application of the present theory is to a pyroelectric plate bonded 
to a fixed foundation, with the upper plane face exposed to sunlight. For this 
the boundary conditions include (i) the prescription of temperature on the 
upper bounding plane, and (ii) the condition of assigned displacement on 
the lower bounding plane. Furthermore, the prescribed boundary values may 
be understood to be functions of a parameter r which depends slowly on 
time: 

T = r(t), \r'(t)\ small. 

Hence, we refer to equations (fT9 l) -( l23l) augmented by these slowly varying 
boundary conditions as a boundary value problem of quasi-statics. 
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4.1 Boundary Control Problem 7.1.3 

4.1.1 Statement of the problem 

The plate V is bounded by the parallel planes x\ = ±h and is coated by an 
infinitesimally thin electrode on the plane x\ = h, so that all its mechanical 
effects may be ignored. We seek solutions of the form 

T = T(x 1 ), <j> = <j>(x 1 ), u i = u i (x 1 ), (24) 

which when substituted in ( TT9|) - (T23T) give 

-/5iTi + c n u ljU = 0, c m u 211 = 0, (25) 

C 4 4M 3> ii + eis^,!! = 0, (26) 

ei5«3,ii + &iT,i - £n0,ii = 0, -%T U + «fi0 )U = 0. (27) 

BCP J. 1.3 : To find the solution of the form (Efy to the field equations 
[25\)-(27\), subject to the ten boundary conditions 

T{h)=T, 0(7i) = 0, h(h)=ti, t 6 (h)=h, t 5 {h)=t 3 , 
w i (-A)=« j (i = 1, 2, 3), -D l (-h)=D, -qi(-h) = q, (28) 

wift T, 0, ti, ?2, ^3) ^1, w 2 , -D and q assigned real constants. 

4.1.2 General solution of BCP 1.1.3 
By Proposition 16.11 with 

c = c 44 , e = e' = eis, = u)i, e = En, k = Ku, k' = K n , (3 = 0, (29) 

K = k/k', A = cu, B = ee' + ce, a = AK~ l B~ l , (30) 
we have from (11001) that the general solution to equations (126"]) - (127]) is 

u 3 (x 1 ) = -c- 1 eKT 1 e ttXl + U 31 x x + U 32 , (31) 

T(xi) = Tie"" 1 + T 2 , (32) 

0(xi) = KT ie aXl + F lXl + F 2 , (33) 

where Ti, T 2 , Ti, F 2 , U31, U32 are arbitrarily chosen smooth functions 
of r. 

Further, (1321) implies Ti = Tiae axi which together with (l23l) 1 yields 

ui(xi) = a" 1 ^! Ti e a:Cl + U u xi + U 12 , (34) 

while ( 1231) 9 gives 

m 2 (xi) = C/21X1 + C/22 , (35) 
where C/ a| g, a, /3 = 1, 2, are arbitrary smooth functions of r. 
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4.1.3 Decomposition of BCPI.l. 3 

We solve BCP J. 1.3 by decomposition into two parts, described below. 
Part 1 of BCP J. 1.3. We first consider the boundary conditions 

T(h)=T, 0(/i) = 0, h(h)=t 3 , 
u 3 {-h) = u 3 , -D 1 (-h) = D, - qi {-h) = q. (36) 

Note that by Eq.s (HHDi, (15T]) - (155]) and (fTUTl) . the 5-t/i boundary condition 
above becomes 

- eit 3) i + £0, i - cuT = -cjT 2 + eFi - eU 31 = D , (37) 

and by ffToT^ the 6-i/i boundary condition above respectively 

kT i - i = -k'Fx = q . (38) 

Now (JSID-dMD satisfy (JMD provided 

T x e ah + T 2 = T , KTx e ah + F x h + F 2 = ^ , (39) 

eiq + c^i = * 3 , -c-'eKT, e~ ah - C/ 31 /i + U 32 = u 3 , (40) 

- ujT 2 - e£/ 3 i + eF l = D 1 -k'Fi = g . (41) 

By solving the above system of equations in the unknowns (Ti, T 2 , Fi, F 2 , C/31, 
we obtain 

r 2 = -££±* - m = e -(T-r,) , (42) 

cu; cxufc V / 

Fi = -Kk- l q, F 2 = 4>-hF 1 - Ke ah Ti , (43) 

1 Fe 
C/31 = -(* 3 - eF x ) , f/32 = u 3 - hU 31 + — -rT x , (44) 
c ce 

so that the solution to the first part of BCP J. 1.3 becomes in particular 

(45) 

= f e »(^)-lV^±^ + ^^e^-^T. (46) 
V / \ ecu CUJK J 



Hence, 



-2ah ,\(m 1 CjD + et 3 . + ^ -A . ~ . 2/i 



-ft) = *( e — - Ij + + A'— 5 j + * + ¥ 5, (47) 
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e -2ah _ 1 



T(-/i) = cD + eh + Kk-^ce + e^q + e~ 2ah T , (4 



and for = D = = q we have 

(f>(-h) = 4> + K{e~ 2ah - ljT, T(-h) = e- 2ah T, (49) 

which yield electric potential and temperature in the plane x\ = —h in terms 
of electric potential <fi and temperature T at x\ — h. 

Part 2 ofBCP 1.1.3. Next, we use ggj, along with ([3T]), (131, fl3g]) 
and ffTTj) to determine the solution that satisfies the four remaining equations 
( 1251) joined to the four remaining boundary conditions 

t^h) = ti, t 6 (h) = t 2 , u^-h) = Hi, u 2 (—h) = u 2 . (50) 

Now by (ED and flM]) we have £i = fi x T x e ah +c n U n - faT and £ 6 = c m U 21 ; 
hence by using ( I32l) -( l33l) the boundary conditions (1501) take the form 



c u U n = (3 X T + t x - faTx e ah , c m U 21 = t 2 , (51) 

- U n h + U 12 = u x - a- 1 c u 1 (3 1 T 1 e~ ah , -U 2X h + C/ 22 = u 2 , (52) 

with T\ given by ( 1421) . On solving ( l5Ti) -( l52l) for U a p and substituting the 
resulting expressions in ( 1341) . ( 1351) we are led to the complete solution of the 
second part of BCP 1.1.3. 

Remark 4.1 We point out that in order to avoid growth as h — > +oo of 
the magnitude of the gradient 

|Vu| + |VT| + |V0| 

of any solution (u, T, 0) to £?CP 1.1.3, we assume 

C44U1 



K(ef 5 + C44SU 



>0. (53) 



In fact, note that by the equalities above, if a < 0, then u^, u 3 , T, <f> — > 00 
as /i — > 00. 

4.1.4 On controllability in BVP J.1.3 

By Eq.s f l4"5"l) . fj4T)]) we can deduce the following control property. 



Remark 4.2 Let — /i < x\ < h. For each choice of Hi, u 2 , TZ3, t±, t 2 , 0, 
given any three quantities in {D, f 3 , 0, T\ , the remaining quantity can be 
choosen to control either T(x\) or 0(xi). 
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4.2 Problem BCP 11.12 



4.2.1 Statement of the problem 

Here V , just as in BCP 1.1.3, is bounded by the parallel planes x\ = 
±h each coated by an electrode, which is infinitesimally thin, so that all 
mechanical effects may be ignored. We seek solutions of the form f[2"4")) which 
when substituted in dTHD"© give again Eq.s (f2oD-(P27j). 

BCP I LI. 3 : To find the solution of the form (24\ ) to the field equations 



[2b°\) -(2l\), which satisfies the ten boundary conditions 

T(h) = T, 4>{h) = 4>, t l (h)=t l , t 6 (h)=t 2 , t 5 (h)=t 3 , 
u i (-h)=u i (2 = 1, 2, 3), 0(-/i) = 2 , -qi(-h) = q, (54) 

with T, <f>, ti, t 2 , t 3 , Mi, u 2 , U3, 4>2 an d q assigned real constants. 



4.2.2 General solution of BCP 77.1.3 

Insertion of ( [241 into the equilibrium field equations (fT9l) -( f23l) gives Eq.s 
(f25l) - ([271) . whose general solution is expressed, as before, by Eq.s ([3Tl) -( l33l . 
where 

T u 7 2 , Fa, F 2 , U al , U a2 (a = 1, 2, 3) 
are arbitrary smooth functions of r. 



4.2.3 Decomposition of BCP 77.1.3 

We solve BCP 77.1.3 by separating it into two parts, described below. 

Part 1 of BCP 77.1.3. We first determine the arbitrary constants in 
the general solution (T3"ll - (!3"3l) so that the boundary conditions 

T(h)=T, (f>(h) = 4, t 5 (h)=t 3 , 
u 3 (-h)=u 3 , (j){-h) = 2 , - qi {-h) = q (55) 

are satisfied. Note that by ({TdT) -?. (!3"2"1) and (!3"3l) the last boundary condition 
becomes 

Kf x KT 1 - K^<p, 1 = -kR- 1 F 1 = q . (56) 

We have 

T x e ah + T 2 = T , KTi e ah + F 1 h + F 2 = , (57) 

eFt + cU 31 = t 3 , -c-'eKT, e~ ah - U 31 h + U 32 = u 3 , (58) 
KT\ e- ah - F x h + F 2 = 2 , -kR- 1 F 1 = q, (59) 



11 



and by solving the above system of equations for the unknowns (T\, T 2 , 7\, F 2 , U 3 
we find expressions for the latter in terms of the boundary data. In particular, 
we have 



2h 



k(e ah - e~ ah ) 



q + 



k' f- -r 



2h 



T 2 = T- e ah T x , 



Q 



-^q- Ke ah Ti 



(60) 



(61) 



Hence, by 
data are 



*i, 



T( Xl ) 



k' k 
the expressions of T and <fi in terms of the boundary 

ah 

, c ~ — c 

2 , 



™ 2A r 

~k 



k! 



2h 



~ q+ 2h 



"-2h^~ 



^ah g— ah 



from which we deduce 

T(-h) = T -2hk- 1 q-k'k~ 1 ( 



(62) 
(63) 

(64) 



Part 2 of BCP 7J.1.3. The remaining two equations, together with the 
appropriate boundary conditions, exactly coincide with BCP 7.1.3. Hence, 
we can proceed as described at the end of Subsection 14.1.31 



Remark 4.3 In order to avoid growth as h 
solution, we again assume (E3J). 



-oo of the magnitude of the 



4.2.4 On controllability of temperature 

By Eq. ( 1621) we can deduce the following control property. 

Remark 4.4 Let —h < x\ < h. For each choice of u\, U2, uz, ti, t 2 , ts, 
given any three quantities from {T, q, <fi, <fi 2 } , the remaining quantity can 
be choosen to control T{x\). In particular, if T and q are assigned, then 
T(xi) is controllable by <p 2 — 0. 



4.3 BCP 7.3.3, plate perpendicular to the polarization 
direction 

Now consider a plate occupied by the same above material but having the 
polarization direction X3 perpendicular to the plane of the plate. The plate 
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is coated by an infinitesimally thin electrode on the plane x 3 = h, so that all 
its mechanical effects may be ignored. Solutions of the form 

T = T(x 3 ) , (j) = <f>(x a ) , Ui = Ui(x 3 ) , (65) 

satisfy (fl9>(p3|) provided 

C44« 1 33 = , C44U2 33 = , (66) 

C 33%,33 - P3T3 + e 3 30,33 = 0, (67) 
e 33M 3 ,33 + ^3 T ,3 - ^330, 33 = 0, -K33T33 + ^3^33 = 0, (68) 

which are a system included in the general case considered in the Appendix. 

BCP 1.3. 3 : To find the solution of the form (E3J] to the field equations 
(E2P-(E2j) ; subject to the ten boundary conditions 

T(h)=T } <Kh) = $, t 3 (h)=t u t A (h) = t 2l t 5 (h)=t 3} 
Ui (-h)=Ui (z = l,2,3), -D 3 (-h) = D, - q3 (-h) = q. (69) 

4.3.1 General solution of BCP 1.3.3 
In particular, on setting in fl99|) 

c = c 33 , e = e' = e 33 , (3 = (3 3 , u = Cb 3 , e = e 33 , k = k 33 , k' = «f 3 , (70) 

we obtain Eq.s fl67|) - fl68|) . Then by Proposition 16.11 the general solution to 

(EZD-dnSD is 

iti(a; 3 ) = U n x 3 + U n , u 2 {x 3 ) = U 2 ix 3 + U 22 , (71) 

u 3 (x 3 ) = a- 1 c- 1 VT 1 e M,s + f/31^3 + C/32 , (72) 

T(x 3 ) = T ie axs + T 2 , (73) 

0(x 3 ) = KT ie aa;3 + F lX3 + F 2 , (74) 

where 7\, T 2 , Fi, F 2 , U a i, U a2 (a = 1, 2, 3) are arbitrary smooth functions 
of r and the notation (1701) and fllOlH is used. Moreover, the constant a, 
given by 

/^3e 3 3 + C 33 ^3 ff7t .s 

a '■= 777 9 r , (75) 

K(elz + c 33 e 33 )' 

is supposed positive. 
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4.3.2 Decomposition of BCP 7.3.3 



We solve BCP 1.3. 3 by decomposing it into two parts. Part 1 of BCP 
1.3.3. First, we note that Eq.s (Tl8>, (I72l)-(I71) and fTTTTO imply 

D 3 (x 3 ) = wTi e aa;3 + uT 2 + eU 31 - eF 1 , (76) 

and consequently, the solutions (172l-(174l meet the boundary conditions ( 1361) 
when 

T x e ah + T 2 = T , KT X e ah + F x h + F 2 = , (77) 
- (3T 2 + eF a + cC/ 3 i = \ , a-^-VTi e" a/l - U 3l h + U 32 = u 3 , 
- ujT 2 - eU 31 +eF 1 = D, -kK~ 1 F 1 = q . 



By solving the above system of equations in the unknowns (Ti, T 2 , Fx, F 2 
we obtain 



etx + cD + Kk- 1 (e 2 + ce)q 



(/3e + cu)e ah 
Q 



Tx 



-ah 



(T - T 2 



Fx = - 



- hFx - Ke ah Tx , 



U 31 = n _ , , ^32 = «3 + ftC/31 - — ^T T 1 > 



ace 



/3e + cu; 

which on substitution in the general solution (fT3l - (174l yields the solution to 
the first part of BCP 7.1.3. In particular, we have 



(x 3 ) = KTxie 



^(1X3 



all 



k'- L q(x 3 - h) + 



T(x 3 ) = T 2 (l - e a{x3 - h) ^j + Te< X3 ~V 
and thus, by (1H01 too, we have 



-h) = K 



n-h) 



T + 



etx + cD + fc /_1 (e 2 + ce) q 
(J3e + cuj)e ah 



e -«_i) + », + ; 



eti + aD + /c' _1 (e 2 + ce) q 



(f3e + cu)e ah 
For = tx = D = 5 we thus obtain 



-2a/i 



1) + e - 2ah 'T 



-h) = K [e~ lah - 1 T + 0, r(-/i) 



-2ah 



T 



78) 
79) 
C/31, C/32) 



80) 



51) 



82) 



83) 
84) 

85) 
86) 

87) 
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which yields electric potential and temperature in the plane x\ = —h in 
terms of electric potential and temperature <ft > T at X\ = h. 

Part 2 of BCP 1.3. 3. The remaining two equations (l66"j) . subject to the 
remaining four boundary conditions ( |50l) . can be solved exactly as before in 
Subsection 14.1.31 

Remark 4.5 We point out that in order to avoid growth as h — > +oo o/ 
the magnitude of the gradient 

|Vu| + |VT| + |V0|, 

of any solution (u, T, 0) to BCP 1.1.3, we assume 

03^33 + C33CJ3 

a := — — ^ > . (88) 

K(e 2 33 + 033^33) 

4.4 On controllability in BCP 7.3.3 

By Eq.s ( 153*1) . ([51]) we can deduce the following control property 



Remark 4.6 Let —h < x% < h. For each choice of Ui, U2, U3, t 2 , t 3 and 
(f), given any three quantities from \t 1 , D, q, T} , the fourth quantity can be 
choosen to control either T(x 3 ) or <p(xs). 

4.5 BCP 77.3.3, plate perpendicular to the polarization 
direction 

Here V (see BCP J. 3. 3) is bounded by the parallel planes 23 = ±h on 
which are coated two infmitesimally thin electrodes whose mechanical effects 
therefore may be ignored. We seek solutions of the form (T5"""l) which after 
substitution in QB-tpi) give Eq.s (pl)-(j?58l). 

BCP J/.3.3 : To find the solution of the form [65°\) to the field equations 
(1Z^-(E2P, which satisfies the ten boundary conditions 

T(h)=T, 0(/i) = 0, t 3 {h)=t 1 , h(h)=t 2 , t 5 (h)=t 3 , 
Ui (-h)=Ui (z = l,2,3), <j)(-h) = 4 2 , -q 3 (-h) = q. (89) 

4.5.1 General solution of BCP IL3.3 

The equations correspponding to f|66l ) - (|68l) are of the form ( "991) on setting 

c = c 33 , e = e' = e 33 , (3 = /3 3 , u> = u 3 , e = e 33 , k = k 33 , k' — «f 3 , (90) 
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so that by Proposition 16.11 the general solution to Eq.s (J6"6"|) - (J6"8|) is given by 
Eq.s (HED-flHD, where 

71, T 2 , 71, F 2 , U al , U a2 {a = 1, 2, 3) 

are arbitrary smooth functions of r and we adopt the notation defined in 

Eq.s (jzoj), (unu. 



4.5.2 Decomposition of BCP 77.3.3 

We solve BCP 77.3.3 as follows: 

Part 1 of BCP 77.3.3. The general solution (f72]l - (f71|l to Eq.s (157)1 - 
satisfies the six boundary conditions 

T(h) = T, <Kh)=f, t 3 (h)=t u 
u 3 (-h)=u 3 , <f>(—ti) = fa , -qi(-h) = q 

provided that 

Ti e ah + T 2 = T , KT ie ah + Fih + F 2 = 
- (3T 2 + eiq + cZ7 3 i = *i , a^c^VTx e~ ah - U 31 h + U 32 
K 71 e~ a/l - F l h + F 2 = 2 , -kK- 1 F 1 = q 
which can be solved for the unknowns 

(71, T 2 , 71, F 2 , U 3 i, U 32 ) 
to give in particular the expressions 



(91) 

(92) 
u 3 , (93) 
(94) 



2h r_ k' 

1 ~ k(e ah - e~ ah ) V q+ 2h 

and ([6OD2, (EI]). Hence, ([73]) and (J7JD become 



*3 



r(*s) 

2/i 



7^ 2/i 
~k 



q + 



k' /-r -r 



2h 



k' f-r -r 



f/ ~2p' 



e a%3 _ gift 



J e ah - c 



^ + 0+f (^-^3) 



(95) 

(96) 
(97) 



(9* 



and we conclude that 

T(-/i) =T-2hk- 1 q-k'k- 1 (J>-^> 2 ) 

Part 2 0/ BCP 77.3.3. The remaining two equations, subject to the 
appropriate boundary conditions, exactly coincide with the corresponding 
equations of BCP 7.3.3. Hence we can solve them by the method of Subsec- 
tion HXS 

Remark 4.7 In order to avoid growth as h — > +00 in the magnitude of 
the gradient of the solution, we assume 
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4.6 On controllability of temperature 

By Eq. ( l96l) we can deduce the following control property. 

Remark 4.8 Let —h < x% < h. For each choice of t and u, given any 
three quantities from {T, q, (p, (j) 2 }, the remaining quantity can be choosen to 
control T(x%). In particular, when T, q are assigned, T(xs) is controllable 
by ~4>-~4> 2 - 

5 Conclusions and perspectives 

We have shown that, for a piezothermoelastic plate referred to a natural 
configuration, in the presence of a quasi-static incremental temperature given 
on one of its bounding faces, on the other bounding face either the electric 
potential or the temperature can be controlled by certain boundary data. 

An aim of a future investigation could be to examine how these result 
generalize when the initial configuration of the plate is not a natural configu- 
ration, that is, when there is some initial mechanical, thermal and/or electric 
field. 

6 Appendix 

The following elementary result on first order differential equations is used. 

Remark 6.1 If f = f(x) is a scalar function of the real variable x and 
a, b G M, then the general solution of the linear first-order differential equa- 
tion 



is f = ^ye ax + b, where 7 is an arbitrary real constant. 

Proposition 6.1 Let c, e, e', (3, 00, e, k, k' be real scalars. Then the system 
of linear differential equations 



f 



a(f - b) 



cu tXX - f3T x + e'0, 
eu tXX +uT x - £0, 
-kT xx + k'4> t 



XX 



XX 



XX 















(99) 



in the unknown scalar functions 



T = T(x), (f> = 4>(x) 



u = u(x) , 
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of the real variable x, has the general solution 

T{x) = 71 e ax + T 2 
<f>{x) = KT ie ax + F lX + F 2 (100) 
u = ar x c- x VT x e^ + U x x + U 2 

where 

(T u T 2 ,F u F 2 ,U u U 2 )eIR 6 

are arbitrary and 

K := k/k', A:= f3e + clu, B := ee + cs, a := AK^B^ 1 , V := (3 - Kae . 

(101) 

Proof. Equation yields (f> tXX = KT XX , thus 

4>,x = KT x + Fi, F 1 eM, (102) 
and Eq.s o become 

u,xx = c~ l {(3T^ x - e'KT xx ) , u xx = e _1 ( - uT x + eKT xx ) . 

By eliminating u jXX from these two equalities, we obtain the second-order 
equation 

T xx = aT x , (103) 

with a defined in ( 11011) : consequently T x = a(T — T 2 ) where T 2 is an arbi- 
trary constant. By Remark 16. II the latter equation has the general solution 

T(x) = T ie ax + T 2 (71, T 2 ) G M 2 , (104) 

which by substitution in (11021) enables us to conclude that (IIOOI) ? holds. 
Lastly, insertion of the expressions for T and into HlU3j) i yields w 
c~ 1 T 1 V ae ax where V := (3 — Ke'a. Hence by integration we obtain m > 
cr l T x Ve ax + Ui, which yields <$M) 3 - 



, XX 
, X 



7 Acknowledgments 

The author would like to thank Professor R. Knops for his discussion and 
suggestions on the present paper. 



18 



References 



[1] Paul, H.S. and Renganathan, K., 1985. Free vibrations of a pyroelec- 
tric layer of hexagonal (6mm) class. Journal of Acoustical Society of 
American. 78, 395-397. 

[2] Paul, H.S. and Raman, K.G.V., 1991. Vibrations of pyroelectric plates. 
Journal of Acoustical Society of American. 90, 1729-1732. 

[3] Eringen, A.C. and Maugin, G.A., 1990. Electrodynamics of Continua 
I. Springer- Verlag New York Inc. 

[4] Nowacki, W., 1975. Dynamic problems of thermoelasticity. Noordhoff 
Int. Publishing, PWN-Polish Scientific Publishers, Warszawa. 

[5] Yang, J.S., 1998. Nonlinear Equations of Thermoviscoelectroelasticity. 
Mathematics and Mechanics of Solids. 3, 113-124. 

[6] Chen, W.Q., 2000. On the General Solution for Piezothermoelasticity 
for Transverse Isotropy With Application. Journal of Applied Mechan- 
ics. 67, 705-711. 

[7] Tiersten, H.F., 1969. Linear Piezoelectric Plate Vibrations. Plenum 
Press, New York. 



19 



